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Abstract 



We compare the normal-state resistivities p and the critical temperatures T c 
for superconducting d x 2_ y 2 pairing due to antiferromagnetic (AF) spin fluc- 
' tuation exchange in the context of the two phenomenological dynamical spin 

susceptibility models, recently proposed by Millis, Monien, and Pines (MMP) 
and Monthoux and Pines (MP) and, respectively, by Radtke, Ullah, Levin, 
and Norman (RULN), for the cuprate high-T c materials. Assuming compara- 
ble electronic bandwidths and resistiviies in both models, we show that the 
RULN model gives a much lower d-w&ve T c (~ 20K) than the MMP model 
(with T c ~ 100K). We demonstrate that these profound differences in the 
. T c 's arise from fundamental differences in the spectral weight distributions of 

the two model susceptibilities and are not primarily caused by differences in 
the calculational techniques employed by MP and RULN. The MMP model, 
claimed to fit NMR data in YBCO, exhibits substantial amounts of spin fluc- 
tuation spectral weight up to an imposed cut-off of 400meV, whereas, in the 
RULN model, claimed to fit YBCO neutron scattering data, the weight is 
narrowly peaked and effectively cut-off by lOOmeV. Further neutron scatter- 
ing experiments, to explore the spectral weight distribution at all wavevectors 
over a sufficiently large excitation energy range, will thus be of crucial impor- 
tance to resolve the question whether AF spin fluctuation exchange provides 
a viable mechanism to account for high-T c superconductivity. The large high- 
frequency boson spectral weight, needed to generate both a high d-wave T c 
and a low normal-state resistivity, also implies large values, of order unity, for 
the Migdal smallness parameter, thus casting serious doubt on the validity of 
the very Migdal approximation on which high-T c spin fluctuation exchange 
models are based. 

PACS numbers: 74.20.Mn, 74.25.-q, 74.25.Fy, 74.25.Ha 
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I. INTRODUCTION 



Recently, an increasing, but not uncontroverted amount of experimental data has been 
shown to be consistent with a d x 2_ y 2 symmetry of the superconducting order parameter in 
the cuprate high-T c superconductors.0'1 In early studies of the Hubbard model, Scalapino 
and co-workersH proposed that a pairing state of d-wave symmetry can arise due to ex- 
change of antiferromagnetic spin fluctuations.! In the model believed to be most relevant 
to the cuprates, the two-dimensional (2D) Hubbard model, a d x 2_ y 2 pairing instability was 
suggested by Bickers et a/.@ and by Moriya et a/J! Recently, two phenomenological models 
of the dynamic spin susceptibility x(<?? for the cuprate materials have been developed by 
fitting an assumed analytical form to experimental data and then using this information to 
calculate such properties as the resistivity p and superconducting transition temperature T c . 

The first model, due to Millis, Monien, and Pines, hereafter referred to as MMP, was 
obtained from fits to nuclear magnetic resonance (NMR) data on YBa 2 Cu 3 7 (YBCofi. 
This model and modifications thereof were then used to calculate T c in weak coupling)! 
and strong couplingHill approaches. Moreover, the resistivity due to spin fluctuation scat- 
tering was calculated, both withoutS and with0 inclusion of leading-order vertex correc- 
tions. Also calculated was the suppression of T c due to non-magnetic impurities.0 The pri- 
mary conclusion drawn from the MMP model calculationsB~0 was that the NMR data, the 
temperature-dependent normal-state resistivity, the large superconducting transition tem- 
perature T C = 100K, and the T c -suppression due to impurities can be consistently accounted 
for with a single, physically acceptable set of model parameters. 

The second model, due to Radtke, Levin, Ullah, and Norman,@ referred to as RULN 
hereafter, was obtained from fits to neutron scattering data on FI?CoJi! This x{q> u ) was 
then used to calculate T c , using a strong coupling approach@'0, the resistivity, by variational 
Boltzmann equation and Fermi-surface-restricted Kubo formalisms^!, and the T c suppression 
due to impurities." The primary conclusion drawn from these calculations was that it is 
essentially impossible to generate a d-wave T c of 90K (or larger) by AF spin fluctuation 
exchange. In fact, if the effective spin fluctuation coupling constant in the RULN model was 
adjusted to reproduce the observed magnitude of the normal state resistivity, the resulting 
d-wave T c was found not to exceed about 10K. 

We are thus faced with the obvious controversy that the two different model fits to 
the dynamic spin susceptibility of YBCO give strikingly different results for the d-wave T c 
and/or the normal-state transport scattering rates. The MMP model allows to obtain a 
very high T c , consistent with the cuprates, as well as a rather low resistivity, consistent with 
recent single crystal data on YBCO. The RULN model, on the other hand, gives a much 
lower upper limit on T c , even for resistivities much larger than the MMP model. 

This controversy has led to a number of criticisms on both sides. An early argumentl^E 
concerned the use of different calculational approaches, specifically the use of weak-coupling 
vs. strong-coupling and approximate semi-analytical Fermi-surface-restricted vs. exact nu- 
merical full-zone approaches to calculate T c . However, both groups@0@ have now carried 
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out detailed strong-coupling, full-zone T c calculations for both models, using fast Fourier 
transform (FFT) techniquestj which allow to include include the full momentum and fre- 
quency dependence of the Migdal-Eliashberg equations. In these full-zone results, the order- 
of-magnitude discrepancy of the T c values between the two models persists. Clearly, the 
different magnitudes of T c are not primarily due to the use of different calculational ap- 
proaches. 

Another objection00 against the original RULN approach concerned their@ use of varia- 
tional Boltzmann or Fermi-surface-restricted Kubo formalisms for calculating the resistivity 
in the RULN model. These transport formalisms, again, amount essentially to an ap- 
proximate Fermi-surface-restricted analysis of the relevant full-zone Bethe-Salpeter integral 
equation. By contrast, in Ref. [Tl|, the resistivity in the MMP model was obtained directly 
from the numerical full-zone solution of the Bethe-Salpeter equation, without additional 
approximations. In the present paper, we will show that, again, the discrepancy between 
the two models is not primarily caused by the use of two different calculational approaches. 

What we wish to demonstrate here is that the very different conclusions about attainable 
T c 's derived from the two models are a direct consequence of the fundamental differences 
in the assumed dynamic spin susceptibilities x{Qi u )- We re-emphasize0 that, even if one 
uses the same electronic band parameters in both models, the RULN model will still give 
a substantially lower superconducting T c than the MMP model. We will show that this 
difference in T c arises, ultimately, from the difference in the spectral weight distribution of 
the two model spin susceptibilities. We therefore conclude that more experimental work is 
needed, primarily by neutron scattering, to elucidate the frequency dependence of the spin 
susceptibility x(<?> u )- Only a complete understanding of the spectral weight distribution, 
including its full frequency and momentum dependence, will eventually allow us to resolve 
the issue of whether or not a spin fluctuation pairing mechanism may constitute a viable 
model of high temperature superconductivity. 

On the theoretical side, our analysis of the two different models shows that, in order 
to produce a large T c and, at the same time, a sufficiently low normal-state resistivity, the 
spin fluctuation models, just like any other boson exchange model, must exhibit substantial 
amounts of spectral weight at a fairly large boson energy scale Q, as well as a substantial 
overall Eliashberg coupling parameter A, assuming typical bandwidths and Fermi energies 
€f in the cuprates. For such large boson energy scales and coupling strengths, the Migdal 
smallness parameter s = XQ/tF is not really small anymore, compared to unity. Thus, 
Migdal's "theorem", the very foundation of current spin fluctuation exchange models, tends 
to be on rather shaky grounds. 

Another important theoretical consequence of having such a large boson energy scale, 
Q/ep ~0.1 — 0.2, is that the conventional mechanism of evading Coulomb interactions by 
retardation becomes largely inoperative in these types of high-energy boson exchange models. 
Thus, even though the local, Hubbard-type on-site Coulomb correlations are completely 
projected away, simply due to the non-s-wave order parameter symmetry, the extended 
(e.g. nearest neighbor) part of the Coulomb interaction potential can be quite effective 
in suppressing the spin-fluctuation-mediated d-wave Current spin fluctuation models, 
including the MMP model,H" have so far entirely ignored the extended part of the Coulomb 
interaction potential and may thus be severely overestimating the actual c?-wave T c values 
which can be ultimately achieved via a spin fluctuation exchange mechanism. 
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In passing, we should also comment briefly on the on-going controversy concerning im- 
purity effects. Monthoux and PinesS claimed that the suppression of T c due to pairbreaking 
by non-magnetic impurities is much weaker than would be indicated by Abrikosov-Gor'kovEi 
(AG) theory. Radtke et a/.0 on the other hand claimed that the suppression of T c due to 
non-magnetic impurities in both models is roughly in agreement with a strong-coupling ver- 
sion of AG theory proposed by Millis et a/H At this point, we can only emphasize, againjl! 
that the two models are indeed quite consistent in that aspect. The differences claimed to 
exist by Monthoux and Pines0 are entirely a matter of different calculational approaches 
for evaluating the impurity scattering rate from the input impurity potential strength, using 
either the i-matrixO or leading-order Born@0 approximation. If the calculated impurity- 
induced T c -suppression is plotted as a function of the calculated impurity scattering rate 
(rather than as a function of the input impurity potential strength), then both approaches, 
Born and t-matrix, give quite similar results. That is, for a given magnitude of impu- 
rity scattering rate, one obtains essentially the same T c -suppression, regardless of how that 
scattering rate was calculated from the input impurity potential. The important point to 
remember here is that it is the impurity scattering rate, not the input potential strength, 
which is measured experimentally via the impurity- induced "residual" resistivity. 

The remainder of the paper is organized as follows: 

In Section II, we introduce the two AF susceptibility models, proposed by MMP0 and 
RULN0, and summarize the different approaches used by Monthoux and PinestMiil (referred 
to as MP, hereafter) and by RULnEI in determining the electronic band parameters and in 
calculating transport properties, for their respective models. 

In Section III, we discuss in some detail the underlying physical assumptions on which, 
implicitly, the different approaches for estimating the band parameters and for calculating 
transport properties are based. We then give a direct comparison of the different transport 
formalisms, for the case of the MMP model, and discuss the limits of applicability of both 
the Fermi-surface-restricted Boltzmann and the full-zone Bethe-Salpeter formalisms, and of 
the underlying Migdal approximation. 

In Section IV, we discuss the origin of the differences in the ci-wave pairing T c 's extracted 
from the two models. We show that the substantial differences in the T c 's of the two models 
are, ultimately, not due to the differences in the assumed bandwidth parameters or in the 
transport formalism used to estimate the coupling parameter. By means of a McMillan- 
Allen-Dynes analysis, we demonstrate that the primary cause for the differences in T c 
lies in the different spectral weight distributions of the two model spin susceptibilities. 

Concluding remarks are presented in Section V. 



II. MODEL SUSCEPTIBILITIES, BAND AND COUPLING PARAMETERS 

We start by describing the two models. The MMP model is a fit designed to describe 
NMR data in YBCO, based on a susceptibility spectral function of the analytical form 
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Im XMM p(g> + .0+) = - * Q ^[ . . - QR-M) (1) 

(1 + - Q| 2 ) 2 + {oj/uj sf Y 

where Q = (n/a, n/a) and q = (q x ,q y ) is in the 1st quadrant of the two-dimensional (2D) 
square lattice Brillouin zone and a frequency cut-off taken as fL, = 400meV. In the most 
recent version of the model, proposed by Monthoux and Pinegiil, hereafter referred to as 
MP-II, the other parameter values are, as given in Table II of that papei0 for a 25% hole 
doping concentration, 

uj sf = UmeV, £/a = 2.3, XQ = 75eV-\ (2) 

where a (=3.8A) denotes the 2D CuOi square lattice constant. Note, that in MP- 110, the 
resistivity as a function of temperature T was actually calculated by using a T-dependent 
^sfi XQi £ as defined in Eqs. (37), (27) and (5b) of that paperj^l respectively, that is, in the 
notation of MP-I10 

u sf (T) = 9.5meV + 4.75 x l(r 2 (meV/K) x T 
Xo^sf 



X Q (T) 



7TU s f(T) 

1/2 



with r s/ = 1.3eV, Xo = 2.6eV" 1 and (3 = 32. At T^IOOK, these T-dependent values Eq.(3) 
agree to within a few percent with those of Eq.(2), taken from MP-II's Table II. In the 
following we will use exclusively the T-dependent values of u s f, £/a, and xq from Eq.(3). 

The RULN model susceptibility, designed to fit neutron scattering data in YBCO, is 
given by0 

ImxRC/L7v(g, u + i0 + ) = C 



1 + J (cos(q x a) + cos(q y a)) 

x 5 3(T + 5V - t -M^ c - M) (4) 

1.05cj 2 -60|cj| +900 + 3(T + 5) 2 V c 1 u KJ 

with the temperature T and excitation energy u measured in meV, a cut-off Q c = lOOmeV, 
a constant J = 0.3 and a prefactor C = 0.19eV _1 . 

In either model, the spin-fluctuation-mediated electron-electron interaction potential 
V(q,u>) is then obtained by multiplying x(g, u;) by a coupling constant g 2 , that is 

V{q,uJ) = g x{<l,u) = g / ; 5 

for transfer of momentum q and complex frequency uo. To calculate superconducting T c 's 
and normal state conductivities, using standard diagrammatic and/or transport theory ap- 
proaches, this interaction potential is combined with a 2D tight-binding electron bandstruc- 
ture 

= — 2ti[cos(ak x ) + cos(ak y ) — 2tcos(ak x ) cos(ak y ) — e] (6) 
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with a chemical potential // = 2e x t± and 1st and 2nd neighbor hopping matrix elements t\ 
and t 2 = tx ti, respectively. To fix the ratios t = t 2 jt\ and e = /i/2t 1; both RULnEI and, 
subsequently, MP- used the Fermi surface shapes and Fermi surface volumes found in 
LDA bandstructure calculations!!, which are consistent with angle-resolved photoemission 
experiments^. In the case of Y Ba 2 Cu? J 0^,j, considered by RULN, this gave t = +0.45, 
working in the hole picture, and a fi value which corresponds to a hole doping concentration, 
measured from half-filling, of a; = 0.18 holes per Cw-site, that is, a total hole concentration 
of rih = 1 + x— 1.18 holes per Cii-site. MP-II adopted the same value for t, but assumed a 
somewhat larger hole doping concentration of x = 0.25 holes per Cu-site, to correspond to 
the case of Y Ba 2 Cu^0^o,. However, MP§§ and RULN0 then employed two very different 
approaches to determine the bandwidth %t\ and the spin fluctuation coupling constant g 2 : 
adjusted t\ to match roughly the calculated LDA bandstructure er near the 
Fermi surface, giving a t\ = 250meV for YBCO. They then adjusted their g 2 to give a 
superconducting <i-wave T c of about 100K, as obtained from solutions of the full-zone k- and 
cu-dependent linearized Eliashberg equations with the spin-fluctuation-mediated interaction 
potential V(q,u) from Eqs.(l), (3) and (5). This gave a value of g = 0.64eV. Using the 
same interaction potential, with the same g 2 and t%, they then solved the full- zone k- and 
cu-dependent Bethe-Salpeter equation in the normal-conducting state to calculate the In- 
dependent resistivity due to spin fluctuation scattering. This calculated model resistivity 
was shown to be in reasonable agreement with the experimental data on YBa 2 Cu^OQ^2. 
Assumed input values and calculated results for the MP-II model parameters are summarized 
in the last column of our Table I. 

RULN, on the other hand, tried to base their parameter values t± and g 2 entirely on 
experimental DC and optical conductivity data, rather than T c and LDA bandstructure 
results. Their starting point are Drude model fits to the measured low-frequency optical 
conductivity of Y Ba 2 Cu^O^,^, of the form 

= "if T+Tjv < 7 > 

parametrized by a Drude plasma frequency ui P; d and a Drude relaxation time Having; 
fixed t = t 2 /t\ and the band- filling, RULN determined their t\ by equating the experimental 
Drude plasma frequency w P) _d, with the "bare" (non- interacting band) plasma frequency 

<o = 2N L ^ J ^ 2 ^M(k) (8) 

where 



Mk) = ^-e 5 , (9) 

a± is the c-axis lattice spacing (perpendicular to the Cu0 2 layers), N L is the number of Cu0 2 
layers per 3D unit cell, and the prefactor 2 arises from the electron spin degeneracy. Note here 
that, for fixed t, a, a±_ and band-filling, uj 2 is proportional to N L xti. Setting u> Pt0 equal to 
the experimental u PtD = 0.97(±0.05)eV for YBCO, RULN thus obtained t x = 80(±10)meV, 
as listed in the column 1 of our Table I. However, in this original estimate, RULN had 
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assumed an incorrect value of Nl = 1, corresponding to only one Cu02 layer per unit cell. 
With the correct number of Nl = 2 Cu02 layers per unit cell in YBCO, one obtains a revised 
estimate for t\ which is reduced from the original value by a factor of 2, i.e. t\ = 40meV, 
since uj 2 (xN L y<ti. Our following discussion will be based on these corrected values for Nl 
and t\, shown in column 2 of Table I. 

Having thus fixed t\, RULN then determined their spin fluctuation coupling constant g 2 
by equating the experimentally determined T-dependent Drude relaxation rate 1/ro, from 
Eq.(7), with the calculated "bare" relaxation rate 1/ro, given by 

l/r = 2ttT x At r , = 2ixT x 2 / — a 2 r F(u) <f>{u/T) (10) 

JO UJ 

where <fr(x) = [(x/2)/sinh(a;/2)] 2 in the variational Boltzmannil and <f>(x) = x/sinh(x) in the 
Fermi-surface-restricted Kubo formalismlllll for x = uj/T. The transport spectral function 
a 2 r F(u), as given by Aliened, is 



2 r , ^ I # J ^(e^Sje^Mk) - ^(fcQ) 2 ImV(fc - j> + z0+) ( 



From To and o; P) o, the DC resistivity in this formalism is calculated as 

P=— = 2 ( 12 ) 

Notice for the following that, as calculated from Eqs. (8-12), for fixed T, t, a and band-filling, 
1/ro and A tr ,o scale proportional to g 2 /h, whereas p scales proportional to g 2 /t\. 

By matching 1/ro, calculated from Eqs. (8-11), in the Fermi-surface-restricted Kubo 
formalism, to a conservative upper limit on the then available experimental 1/ro data in 
the cupratesi , taken between about 120K and 350K, RULN thus inferred an upper limit 
for the coupling constant, g 2 = 0.533eV 2 , based on their incorrect t\ estimate of 80meV, as 
shown in the column 1 of Table 1.00 Using the corrected value of t 1 = 40meV, this g 2 value 
also has to be corrected downward, by a factor of 2, to g 2 = 0.267eV 2 , as shown in column 
2 of Table I, in order to maintain the same values of l/r and p as in the original RULN 
paper. Note that, given the band parameters, the foregoing procedure of determining g 2 
is essentially equivalent to adjusting g 2 so as to roughly fit the calculated resistivity from 
Eq.(12) to the (upper limit) experimental resistivity 

p = 1/Rea(0) = — . (13) 

U ptD T D 

From standard transport theoryll'il it is well-known that neither uj Pj q and uj Pj d nor tq and 
td are necessarily the same. RULN's basic approach of simply equating first u Pt0 = u Pt D, and 
then t = td, to extract the bandwidth 8ti, and, respectively, g 2 from the optical conductivity 
data, is thus based on several simplifying assumptions. This will be discussed in more detail 
in Section III of this paper. 

By solving the full-zone k- and cu-dependent linearized Eliashberg equations for the RULN 
model, with g 2 and t\ values extracted from the optical data, one then obtains an upper-limit 
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ci-wave T c , from Fig. 1, of about 11K, as shown in columns 1 and 2 of Table I. An important 
point to notice here is that RULN based their parameter estimates on the then available^! 
transport data on early YBCO samples@ which had much larger in-plane resistivities than 
the more recently reported samples.E3 Consequently, RULN's model parameters give a re- 
sistivity for their mo which is significantly larger than MP-II's p values for the MMP 
model or the recent experimental p data@ on Y Ba<iCuO§$, used by MP-II for comparison 
to theii0 MMP model results. 

Specifically, using RULN's u) p $ = 0.97eV and theirB upper limit for \t r ,o = 1/(27tTto) ~ 
0.52 (approximately constant for T>120K) RULN have a relaxation rate 1/tq = 2TrX tr ^T ~ 
56meV and a resistivity of about 450/if2cm at T = 200K, as shown in columns 1 and 2 of 
Table I for the Fermi-surface-restricted Kubo formalism. By contrast, MP-II obtained about 
6ApQcm at 200K, from their solution of the full Bethe-Salpeter transport equation, consistent 
with the in-plane-averaged resistivity of Y Ba 2 Cu30Q^. From the width (HWHM) of MP- 
II's Recr(c<j)-results at 200K, shown in Fig. 13 of MP-II, we can also infer a rough estimate 
for the calculated Drude relaxation rate 1/td and Drude plasma frequency ou p ,d in the MMP 
modelO, namely l/ro~12meV and, from Eq.(13), ui P; d = / pro) 1 / 2 ~ 1.2eV, as shown in 
the last column of Table 1.0 

Thus, for the respective model parameters used, MP-II and RULN obtain or, respec- 
tively, assume comparable values for the Drude plasma frequency, u PtD ~ leV, consistent 
with the experimental data. However, the Drude relaxation rate and the resistivity of MP- 
II are smaller than those of RULN by factors of about 5-7. And yet, the <i-wave T c of MP-II, 
is larger than that of RULN by a factor of about 9. We are thus called upon to explain 
why T c in the RULN model comes out so much smaller than in the MMP model, even for 
parameter values giving the RULN model a substantially larger resistivity than the MMP 
model. 



III. MASS RENORMALIZATION, TRANSPORT FORMALISMS, AND MIGDAL 

PARAMETERS 

In the present section, we will discuss in some detail the major calculational differences, 
between RULN and MP-II, in their determination of the bandwidth and in their formalism 
for obtaining the resistivity. We will elucidate the underlying assumptions upon which the 
different approaches are based and point out some of the potential short-comings of either 
approach. 

First of all, regarding the bandwidth determination, recall that the RULN estimate for 
the bandwidth 8t\ = 0.32eV, corresponds to a non-interacting plasma frequency u Pi0 which 
equals the observed Drude value for YBCO, u p<D = 0.97eV (see columns 1 and 2, Table I). 
This bandwidth is a factor of 6.25 smaller than the LDA-based 8t\ = 2.0eV of MP-II. The 
latter corresponds, by Eq.(8), to a non-interacting plasma frequency ui^^ =2.42eV, at the 
RULN's assumed band-filling = 1.18 (see column 3, Table I). It is presently not known 
whether the observed optical mass enhancement 
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JD =(u% DA) /^ D ) 2 -6.25 (14) 

is caused primarily by spin fluctuations or by non-spin-fluctuation effects. It seems likely 
that both contributions are present and non-negligible in the cuprates, that is 



J D — ^nsf 



x Z sf (15) 



with Z s f denoting a spin fluctuation and Z ns f a non-spin-fluctuation contribution, respec- 
tively. In equating the observed Drude plasma frequency uj p ,d with their "bare" model 
plasma frequency u> P fi, RULN thus implicitly made the following simplifying assumptions: 

Firstly, they assumed that all the non-spin-fluctuation-induced mass enhancement effects 
can be essentially accounted for by using, as input into the spin fluctuation calculation, an 
effective bandwidth 8ti and an effective spin fluctuation coupling constant g 2 into which all 
non-spin-fluctuation renormalization effects have been absorbed. Thus 

ry _ ( (LDA) I \ 2 

^nsf — \Upfi I ^p,0 J 

Z s f = (u pfi /uJp,D) ■ (16) 

where uj P) q is the effective "bare" plasma frequency, corresponding to RULN's effective t\ via 
Eq.(8). 

Secondly, RULN assumed that the observed Zd is dominated by non-spin-fluctuation 
effects, i.e. arises from local Coulomb (and possibly other, such as electron-phonon) inter- 
actions excluding the effect of spin fluctuations, whereas the spin fluctuation contribution is 
negligible, that is 

Zd — Z ns f 

|^/-1|<1 (17) 

Given these two assumptions, the effective bandwidth and effective coupling constant can 
be inferred directly from the experimental data for uJ Py D an d 1/td, as was done by RULN. 
Note however, that the term "bare," in the context of the RULN approach, should then be 
understood to mean only "bare of renormalizations due to spin fluctuations", and not "bare 
of all residual interactions beyond LDA" . 

Recall here that, in conventional Fermi-surface-restricted transport theory,Eil!l the width 
(HWHM) of the Drude peak, 1/td, and its integrated spectral weight, uj 2 d , are related to 
t and ujpQ as follows: 

^l,D — ^Ifl/Zsf 
r D = r x Z af . (18) 

Thus, 1/td and u 2 D are both renormalized by Z s f in such a way that Z s f cancels out in the 
DC resistivity 

4W 471 (19) 
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as already implied by Eqs.(12) and (13). One can therefore indeed calculate the DC resis- 
tivity p = l/(x(0) directly from u 2 and 1/t , as assumed by RULN in Eq. (12) .El However, 
one may not equate u> Pt0 with the observed oj p ,di nor r with r D , unless the conditions of 
Eq.(17), i.e. \Z s f — 1|<C1 are met, a point on which we elaborate further below. 

As far as spin- and non-spin-fluctuation mass enhancements are concerned, take 
the extreme opposite point of view: They0 implicitly assume that essentially all of the 
observed Drude mass enhancement Zd is generated by the spin fluctuations and that all 
non-spin-fluctuation effects are negligible, that is, in Eq.(15), 



Given that assumption, one should compare the experimental data only to the calculated 
model Drude relaxation rate 1/td and calculated Drude plasma frequency u Pt D, which fully 
include all spin-fluctuation-induced mass renormalization effects. To estimate the values 
of 1/td and uj p .d given in Table I for the MMP model,0 we have thus extracted 1/td 
directly from the Drude peak width (HWHM) of the Bethe-Salpeter solution for the optical 
conductivity <j(u), shown in Fig. 13 of MP-II.0 Using the p value (at T = 200K) from Fig. 
12 of MP-II, we have then obtained uj Pj d from Eq.(13), i.e. by u> Pt D = (47r/pr d) 1 ^ 2 - The 
resulting MMP model value for LU Pt £>, as discussed above, is in rough agreement with the 
experimental data on YBCO. 

Physically, both the RULN and the MP treatment of the auasi-particle mass enhance- 
ment raises some concerns. As far as the RULN approach" is concerned, their model 
assumption of simply absorbing all non-spin-fluctuation effects into an effective t\ and an 
effective g 2 clearly needs to be re-examined and justified within the framework of a mi- 
croscopic theory. Also, for the t\ and g 2 estimates given in column 2 of Table I, the spin 
fluctuation contribution Z s f to the overall mass enhancement Zd may not really be neg- 
ligible. Namely, given the spin-fluctuation-mediated interaction potential, Eq.(5), one can 
get a rough estimate of Z s f from the self-consistently calculated single-particle self-energy 
T>(k,u + i0 + ) due the spin fluctuations as follows: 



where (...)fs denotes a suitable Fermi surface average and Xz is the dimensionless Eliashberg 
coupling parameter 



obtained from the Eliashberg spectral function of the spin-fluctuation-mediated interaction 
potential 



Zd — Z s f 
Z nsf -1\ < 1 



(20) 






(22) 



a 2 z F{u) 



(23) 
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In the RULN model with the parameter values of columns 1 or 2, Table I, and temperatures 
T~200K one finds Az~0.8 — 1.1 and thus Z s f~ 1.8 — 2.1 which is not a negligible factor, 
contrary to Eq.(17)£a 

To infer both t\ and g 2 from the transport and optical data, one should therefore really 
carry out a calculation for the full dynamical conductivity cr(u) in the RULN model, in- 
cluding all spin-fluctuation-induced mass-renormalization effects, and adjust both ti and g 2 , 
simultaneously, such that the calculated 1/td, from the Drude peak width, and the calcu- 
lated u Pj [), from u Pt D = (4tt / pro) 1 ' 2 , Eq.(13), both match the experimental values. Such an 
approach would likely give substantially larger estimates for both t\ and g 2 if one assumes 
as input RULN's original experimental w p ,d an d (upper-limit) p values, as given in columns 
1 and 2 of our Table I. 

However, if one uses the more recent, much lower p data for YBCO^, the original 
RULN procedure would give a much smaller g 2 and \z estimate, with about the same ti, 
as in column 2 of Table I. The resulting Z s f would then indeed be close to unity, i.e. the 
conditions of Eq.(17) would be satisfied. So, based on the more recent, lower p data, the 
original RULN approach@ could actually be justified, at least as far as consistency with the 
underlying assumptions about Z s f is concerned. It is obvious from Fig.l, that the resulting 
T c , based on these lower p values, would be substantially below 11K, due to the smaller g 2 
value. So, if one were to use the RULN model susceptibility together with the most recent 
transport data@ for YBCO as input into the a fully self-consistent calculation of ct(uj), 
including all spin-fluctuation-induced mass-renormalization effects, one would ultimately 
find about the same t±, but a substantially smaller value for g 2 and T c than the original 
RULN estimates stated in columns 1 and 2 of our Table I. 

As far as the MP approach00 is concerned, there is no good justification for assuming 
that non-spin-fluctuation contributions to the mass enhancement are negligible. As recently 
demonstrated explicitly for the case of electron-phonon coupling,0 non-spin-fluctuation con- 
tributions can significantly suppress the (i-wave T c , if \Z ns f — 1| ~(9(1). In the MMP model, 
one is then forced to assume a substantially larger spin fluctuation coupling g 2 than was 
used in MP-II, in order to maintain a T c ~ 100K. This, in turn, will cause a larger model 
resistivity in the normal state, above T c . Thus, if substantial non-spin-fluctuation contri- 
butions to the the mass enhancement are present, the agreement between the MMP model 
results and the experimental transport data becomes somewhat questionable. 

A second major difference between MP00 and RULN0 lies in their calculational ap- 
proaches for obtaining the resistivity. Originally, MP had calculated the optical conductivity 
cr{ui) diagrammatically from the self-consistent single-particle Greens function, in the sim- 
plest single-loop approximation, without vertex corrections.0 Subsequently, in MP-II,0 they 
tried to improve upon that approach by including vertex corrections to the current correla- 
tion function, at the level of a Bethe-Salpeter ladder approximation.!!! Their solution of this 
"full-zone" Bethe-Salpeter equation was done numerically, without further approximations, 
including the full k- and tu-dependence of the current vertex function. 

RULN, on the other hand, relied on the variational Boltzmann and Fermi-surface- 
restricted Kubo formalisms@Jiiii outlined above. It is well-knownilil that the semi-classical 
Boltzmann equationEa can be derived, under certain simplifying conditions, by an approx- 
imate Fermi-surface-restricted analysis of the full-zone Bethe-Salpeter equations. These 
simplifying conditions are indeed well obeyed by conventional wide-band electron-phonon 
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systems.il'!! However, as discussed further below, it is not clear whether these simplifying 
conditions are also satisfied for the spin fluctuation models under consideration here. 

To test the various transport formalisms, we have calculated the resistivity for the case of 
the MMP model, using both of the Fermi-surface-restricted approaches^ and compared 
to the results obtained by MP-II0 from numerical solutions of the full- zone Bethe-Salpeter 
equation, as shown in their!! Fig. 12, with the T-dependent spin fluctuation parameters 
given by Eq.(3). The results of this comparison are shown in our Fig. 2. As suggested in 
MP-II, we find that both the variational Boltzmann and the Fermi-surface-restricted Kubo 
formalism, as used by RULN, do indeed overestimate the resistivity, compared to the Bethe- 
Salpeter results. On the other hand, the overestimate, being about a factor of 1.7 for the 
Fermi-surface-restricted Kubo formalism at 200 K, is by no means as large as was implied 
m MP-II.0 For the variational Boltzmann formalism, the overestimate is somewhat larger, 
about a factor of 2.1 relative to the Bethe-Salpeter results at T = 200K. Note here, that the 
g 2 estimate obtained by RULN for their model!l was based on resistivity calculations in the 
Fermi-surface-restricted Kubo formalism. 

Several causes could, potentially, contribute to the differences between the full-zone 
Bethe-Salpeter and the Fermi-surface-restricted variational Boltzmann results shown in Fig. 
2. First of all, there is a numerical issue, having to do with the method used in MP-II 
to calculate the optical conductivity ct{uj) in the real-frequency domain. In MP-II, the 
numerical solution of the Bethe-Salpeter equation and the subsequent calculation of the 
current-current correlation function C(u) = iuoa{uj) were actually carried out in the Matsub- 
ara imaginary-frequency domain. The Vidberg-Serene@ Pade-approximant technique was 
then used to analytically continue the imaginary-frequency data C(iuj m ), from imaginary 
frequencies iu m = 27rmiT into the real- frequency correlation function C(u + i0 + ), defined on 
the real-u; axis. This analytical continuation procedure is mathematically equivalent to an 
inverse Laplace transform and well-known to be numerically highly unstable. That is, even 
very small numerical errors in the imaginary-frequency input data C(iu m ), for example due 
to finite cut-offs in the Matsubara frequency summations, can produce large errors in the 
real-frequency output data C(u + i0 + ). This problem tends to become particularly severe 
at higher temperatures when the smallest non-zero Matsubara frequency off the real axis, 
iui, becomes comparable or larger in magnitude than the characteristic frequency scale one 
is trying to resolve on the real-u; axis. For example, in the MP-II calculation at T = 200K, 
u>i = 2ttT = 108meV which is more than 4 times larger than the width of the real-frequency 
spectral feature, the Drude peak width 2/r D = 24meV, to be resolved. It is not clear what 
level of numerical accuracy was actually achieved in the analytical continuation of MP-II.0 

Another potential source for discrepancies between variational Boltzmann and full- zone 
Bethe-Salpeter approach lies in the two-step approximation which is being used to derive the 
former from the latter.!!'!! The first approximation step consists of restricting the analysis 
of the full Bethe-Salpeter ladder equation to the Fermi surface.@ll This step reduces the 
Bethe-Salpeter equation to the linearized semi-classical Boltzmann transport equation. The 
second step is then to employ a simple variational ansatz to obtain an approximate analytical 
solution to the Boltzmann equation, which leads to our Eqs. (8-12). 

The error from the second approximation step, the variational ansatz, may become ap- 
preciable when the temperature falls below the characteristic boson energy scale Q. For 
example, in the case of acoustic phonon exchange, the error is found to cause overestimates 
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in p as large as 30% at temperatures of about | of the Debye phonon energy!!. In the case 
of the MMP model, the relevant boson energy scale Q is of the order of 100 — 150meV, as 
estimated e.g. by the spin fluctuation spectral moments (uj p ) 1 J p listed in the last column of 
Table I and defined below. This is about a factor of 6 — 9 larger than T at T = 200K and 
could therefore cause errors due the variational ansatz to be of similar (~30%) magnitude. 
However, in the case of the RULN model, the respective frequency moments (see Table I 
columns 1-4) are about 2.5 — 3 times smaller, thus substantially reducing the error due to 
the variational ansatz for that model.il 

The first approximation step, reducing the Bethe-Salpeter to the linearized semi-classical 
Boltzmann equation, requires that the dimensionless Migdal parameter 

s = Xz- (24) 

be small compared to unity. Here, Xz is the Eliashberg parameter, from Eq.(22), ep is the 
Fermi energy, measured from the nearest band edge, and Q is again the characteristic boson 
energy scale. A reasonable order of magnitude for Q can be estimated e.g. by taking a 
low-order frequency moment of the Eliashberg spectral function a z F(u), Eq.(23), or of the 
transport spectral function a^ r F(u), Eq.(ll), that is, e.g. 

n ~ (^)f (25) 

where 

(^)z = r -a%F{u)u* I r -alFico) (26) 

JO LO / Jo id 

The results for (oo p ) z p with p= 1 and 2 at T = 200K are given for both models in Table I. 
The moments {lj p )H p of the transport spectral function af r F, from 

= r -alF{u)u p I r -alF(u\ (27) 

Jo U! J Jo 

are not listed in Table I, but are, in the RULN case, identical to the analogous moments of 
a 2 z F and, in the MMP case, within 1% of the analogous o? z F moments. 

In conventional wide-band metals, with transport scattering dominated by phonons, we 
have Az~ 1, f2 of the order of the Debye energy and hence s of order 10 -2 . The Fermi-surface- 
restricted Boltzmann theory is thus expected to work quite well in such systems. However, 
in the MMP model, both a z F(u), from Eq.(23), and a% r F(u), from Eq.(ll), exhibit a wide 
peak, centered around u m ~90 — 120meV (depending on T), with substantial spectral weight 
extending up to the cut-off energy of Q c = 400meV, as shown for T = 200K in Figs. 3a and 
3b. Withe F ~4t 1 ~1.0eV (near ^-filling), A z ~1.6 and fi~100-150meV (fromEq.(25) with 
p~l — 2), the Migdal parameter s at 200K is about 0.16 — 0.24 here. Noticeable quantitative 
discrepancies between the full-zone Bethe-Salpeter and Fermi-surface-restricted Boltzmann 
theory can therefore be expected. 

For the case of the RULN model, both a^ r F(uj) and a z F(to), as shown in Figs. 3a and 
3b, have a sharp peak around o; m ~30 — 50meV and the spectral weight above u m decreases 
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more rapidly and is cut off altogether at fi c = lOOmeV. This is reflected in a substantially 
smaller boson energy f2~30 — 50meV, from Eq.(25) with p~ 1 — 2. For the band and coupling 
parameters from column 2 of Table I, giving an Eliashberg parameter A^~0.8 and a Fermi 
energy e F ~ Ati = 0.16eV, the RULN model Migdal parameter at 200K, s ~ 0.15 — 0.25, is 
comparable to that of the MMP model. 

However, if one chooses t\ and g 2 in the RULN model so as to give the same bandwidth 
and the same resistivity as in MP-II, one finds the parameter estimates given in column 3 
of Table I and discussed further below. For this parameter set, with ep ~4ti = l.OeV and 
Xz ~ 1.3, the Migdal parameter s at 200K is only of the order 0.04 — 0.07, i.e. noticably 
smaller than in the MMP model. For this last RULN parameter set (column 3 of Table 
I), we therefore expect the Fermi-surface-restricted Boltzmann theory to provide a better 
approximation to a full-zone Bethe-Salpeter calculation of the resistivity. 

The overall magnitude of the Migdal s parameter in the MMP model raises serious 
concerns about the validity of the entire spin fluctuation exchange approach0'0 itself. It 
is well-known,ilil0 that the validity of the Migdal-Eliashberg approximation rests on the 
smallness of that same parameter s, Eq.(24), which is also required to ensure the validity of 
the Fermi surface-restricted Boltzmann approach. In other words, the Migdal parameter s 
which controls the relative magnitude of the discrepancy between full-zone Bethe-Salpeter 
and Boltzmann theory transport results, Fig. 2, also controls the relative magnitude of 
those higher-order non-Migdal vertex corrections which have been discarded, both in MP-II 
and in RULN, in the single-particle self-energy, in the particle-particle interaction kernel 
of the linearized Eliashberg T c -equations, and in the particle-hole interaction kernel of the 
Bethe-Salpeter equation. One should therefore not be mislead into believing that numerical 
solutions of the full-zone Bethe-Salpeter equations necessarily constitute a "better" trans- 
port theory than the semi- classical Fermi-surface-restricted Boltzmann approach.iHl Either 
approach works well when s is sufficiently small compared to unity. And either approach 
can be expected to give, at best, only qualitatively correct results when the s parameter 
attains magnitudes as large as estimated above. Concerns regarding the applicability of 
the Migdal approximation in spin fluctuation exchange models have already been raised in 
several earlier and recent studies. 



IV. THE EFFECT OF SPECTRAL WEIGHT DISTRIBUTION ON T c 

The differences in the assumed band parameters and in the transport formalisms of 
RULN and MP-II certainly contribute to differences in their respective <i-wave T c values. 
For the purpose of comparing the spin-fluctuation-mediated pairing in the two models, it is 
therefore of interest to treat the two models on an equal footing, as far as band parameters 
and transport formalism are concerned. In order to demonstrate that the differences in the 
T c values of the two models is not primarily caused by the differences in band parameters 
or transport formalism, we have thus carried out a T c calculation for the RULN model with 
the following modified parameter set, shown in column 3 of Table I: 
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Firstly, instead of ti_ = 40meV, we use the same LDA-based value t± = 250meV which 
was assumed by MP^'tll Secondly with that t\ value fixed, we adjust the RULN coupling 
constant g 2 so that the resistivity, calculated in both models at 200K in the Fermi-surface- 
restricted Kubo formalism Jill! is the same, namely p= 109/if2cm. Recall that in the Fermi- 
surface-restricted Kubo formalism, Eqs.(8-12), the resistivity scales proportional to g 2 /t\, 
for fixed t = t 2 /^i, band-filling and T. Consequently, the modified g 2 value of 2.53eV 2 , thus 
obtained, in column 3, is about a factor of 9 larger than the g 2 estimate in column 2, in spite 
of the fact that the underlying p value is almost a factor of 4 smaller than that assumed in 
column 2. 

The important point to notice now is that the resulting T c = 18. 4K for the modified 
RULN parameter values, in column 3, is increased only modestly over the original RULN 
estimate of 11. 2K. Thus, in comparing the column 3 results for the RULN model to the 
MP-II results, in the last column of Table I, we conclude that the RULN T c value is about 
a factor of 5 lower than the MP-II value, even if we assume the same model bandwidth, 
similar band-filling, and coupling constants g 2 giving similar values for the DC resistivity, 
calculated in both models, by the same transport formalism. 

In fact, as noted earlierOil and clearly seen in Fig. 1, T c as a function of g 2 , in the RULN 
model, does not increase indefinitely with increasing g 2 , but rather seems to approach a finite 
saturation value for g 2 — > oo. This saturation T c does not exceed values of about 50K for the 
band parameters considered here. Using transport data and/or LDA bandstructure results 
to obtain actual estimates for g 2 and t\ only serves to make this fundamental upper limit 
on T c in the RULN model more stringent. Thus, given the RULN spin fluctuation spectrum 
of Eq.(4), with the spectral parameter values stated, it is simply not possible to generate, 
within a conventional Eliashberg theory, a <i-wave T c of, say, 90K or larger, regardless of the 
choice of band parameters and coupling strength g 2 . 

Why then does the the RULN spin susceptibility give such a low T c , compared to the 
MMP model, even when its assumed bandwidth and resistivity is comparable to that of 
MP-II ? To answer this, we define a pairing spectral function as follows 



a 2 J 70 / £&8(e e )5{e2) V (k) V (V)Tm Vi k //. ,s /()- ) 



<f{u) = — ;;;\; A : : ( 2s , 

where 



rj(k) = cos(k x a) — cos(k y a) (29) 

is the basis function for nearest neighbor d x 2_ y 2 pairing. Note here that the sign in Eq.(28) 
is defined such that a positive o? d F corresponds to attraction in the (i-wave pairing channel. 
As shown in Fig. 3c, the spectral weight distribution of a 2 d F is quite similar to that of a 2 z F 
and af r F in either model. Note here that, again, in the MMP model a d F has a very broad 
peak, at an tu m ~35 — 55meV (depending on T), with substantial spectral weight extending 
out to the 400meV cut-off. And, again, in the RULN model the spectral weight is much 
more narrowly peaked, at similar peak position u m ~ 30 — 50meV, and falls off rapidly well 
before the lOOmeV cut-off. We emphasize that this latter spectral weight distribution is 
claimed to be in accordB'i with the currently available neutron scattering data on YBCO. 
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Recall now, from Eqs.(8-ll), that the transport relaxation rate 1/tq, and hence the 
resistivity p, is most sensitive to the low-frequency part of af r F, up to w ~ T. Thus, with 
their u p values being about equal, the two models will give roughly the same resistivity 
if their respective coupling constants g 2 have been adjusted so that their ot 2 r F are roughly 
the same, in absolute magnitude, at low frequencies uj ~ T. However, given such coupling 
constant values, the MMP model will then give a much larger superconducting T c than the 
RULN model, because of its substantially larger spectral weight in a 2 d F at high frequencies. 
The basic physical reason for this is that, unlike the normal state DC resistivity p, the 
superconducting transition temperature T c is quite sensitive to the high-frequency part of 
the boson spectrum. 

To analyze this more quantitatively, we consider the moments and the dimensionless 
Eliashberg coupling parameters associated with 

(u?) d = r -a 2 d F{u)u p I r d -a 2 d F{u) (30) 

JO OJ / JO OJ 



and 



2 r d ^a 2 d F{u). (31) 

JO id 



Notice that all frequency moments defined here and above ((oj p ) 1 J p , (u p )}l P and (oj p ) d ) are 
independent of the coupling strength g 2 and bandwidth 8tx, assuming fixed T, band shape 
t = tijt\ and band-filling. All the A parameters (X z , Xt r ,o, Ad), on the other hand, scale 
proportional to g 2 /ti, at fixed T, t and filling. 

In approximate T c -formulas of the McMillan- Allen-Dynes varietyEHil, T c is given in terms 
of Xz, Ad and (oo p } 1 J p by an expression of the form 

T c = (u; p ) d /p F(X z , X d ) (32) 

where p can be weakly A-dependent with + < p < 2. F is monotonically increasing as a 
function of A^ and monotonically decreasing as a function of Xz- In the weak-coupling limit, 
Ad, Az<Cl, F varies roughly exponentially with Xz and l/Xd, that is 

F(X Z , Xd) = const x exp(-l/A*). (33) 

where 

* ^ TTT Z - < 34 > 

Note however that, for the coupling strengths discussed here, in both models Xd, A#~0(1), 
as shown in Table I. So both the MMP and the RULN model are really in a regime of 
intermediate, not weak, coupling. In this intermediate-coupling regime, the variation of T c 
with Xz and l/Xd is typically less rapid than implied by Eq.(33).c°lcll 

The values of (oj)d, (uj 2 ) 1 / 2 , X d and X d , evaluated at the respective d-wave T c 's, are shown 
in Table I. In comparing the RULN results, from e.g. column 3, to the MMP results, in the 
last column of Table I, there are three key points to be noted: 
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The first point is that the low-order frequency moments of a 2 d F(uj) in the MMP model 
are again significantly larger than those of the RULN model, by a factor of about 2-3. The 
second point is that also the pairing A-parameter, A^, in the MMP model is about 2.2 times 
larger than the RULN value in column 3, despite the fact that both models have about 
the same \t r ,o at 200K. The third point is that, per A^, the MMP model has a noticeably 
smaller mass-enhancement A-parameter, Xz, than the RULN model, that is, the ratio Xz/Xd 
is about 0.88 in the MMP model, compared to 1.74 in the RULN model. The last two points 
imply that X d in the RULN model is noticeablv smaller than in the MMP model. 

In the McMillan-Allen-Dynes formalism,c3EJJ each of the foregoing three differences be- 
tween the two models will tend to give the MMP model a larger T c than the RULN model. 
The first two differences, in (u p ) 1 J p and A^, are an immediate consequence of the substantial 
additional high-frequency spectral weight in the MMP model, for lo ~ 50meV, as shown in 
Fig. 3c. The third difference, in the A^/A^-ratio, is more subtle and arises from differences 
in the detailed momentum dependences of the two model susceptibilities. 

To demonstrate more explicitly that the difference in the high-frequency spectral weight 
is indeed the primary cause for the differences in the T c 's of the two models, we now consider 
what happens if we take, for example, the basic functional form of the RULN model, but 
allow it to have a substantially larger boson frequency scale than originally imposed by 
RULN on the basis of the neutron scattering data.@ll! That is, we modify the RULN spin 
susceptibility model by re-scaling its energy spectrum, leading to a modified susceptibility 
function 



Xrulnj (q, uj) = Xruln {q, w/f) (35) 

with a boson energy scaling factor / and Xruln from Eq.(4). A re-scaled MMP model 
susceptibility Xmmpj could be defined analogously with the MMP suceptibility Xmmp from 
Eqs. (1) and (5). Notice here that all frequency moments, ((jj p ) 1 J p , {oo p )\i P and (uj p )]/ p } scale 
linearly with /, whereas Xz and A^ are independent of /, assuming T and all other model 
parameters are fixed. From full-zone Eliashberg calculations, we find that, upon boson 
energy re-scaling at fixed A parameters, T c varies indeed roughly linearly with (uj p ) 1 J p in 
either model, consistent with the basic McMillan- Allen-Dynes approach, Eq.(32). Thus, by 
increasing the boson energy scale, with f> 1, at fixed g 2 , we can raise the rf-wave T c in the 
RULN model to values which are as large as those in the MMP model, or larger, while at the 
same time lowering the model's normal state resistivity, above T c , from Eqs. (8-12). Likewise, 
by lowering the overall boson energy scale, in either model, with / < 1, we would decrease 
T c and increase the resistivity p. Notice that At ri o and p, from Eqs. (8-12), decrease with 
increasing boson energy scale, since increasing /, at fixed T, actually reduces the relevant 
spectral weight of a 2 r F(uj) at low frequencies u;~T. 

Suppose now we adjust both g 2 and / in the re-scaled RULN model, Eq.(35), such that 
both its resistivity p at T = 200K (from the the Fermi-surface-restricted Kubo formalism, 
Eqs. (8-12), say) and its ci-wave T c (from full- zone Eliashberg solutions) match those of 
the MMP model, as given in the last column of Table I. The results for this last RULN 
parameter set are shown in column 4 of Table I. The required value for the scale factor 
is about / = 2.70, implying that, relative to column 3, all the various RULN frequency 
moments are increased by that factor. Also, the coupling constant g 2 and, hence all A 
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values, are substantially increased, by more than a factor of 2, relative to those of column 
3. This increase in g 2 is necessary in order to maintain a constant resistivity, p— 109/if2cm, 
i.e. in order to compensate for the decreasing effect on p caused by raising /. 

In comparing the RULN results from column 4 to the MMP results in the last column 
of Table I, the crucial point to notice is that now both models, having the same T c and p, 
also have roughly, to within 20 — 30%, the same frequency moments. The fact that now the 
(re-scaled) RULN model acquires somewhat (~ 20 — 30%) larger frequency moments than 
the MMP model can be rationalized by noting that the smaller in the RULN model has 
to be compensated for by a larger boson energy scale, in order to get the same T c as in the 
MMP model. 

From the foregoing results, it is clear though that the spin fluctuation spectral weight 
distribution is the central factor determining T c in either model and that the differences in T c 
are not primarily caused by differences in bandstructure parameters or transport formalism. 



V. SUMMARY 

We can see now that the problem of whether high-T c superconductivity is explained 
by a spin fluctuation model is closely tied to the question about which type of dynamic 
susceptibility is more realistic, MMP or RULN. One important feature of the MMP model 
susceptibility is that its behavior is similar to an effective RPA form which is known to fit 
Monte Carlo data on the single band Hubbard model.0 That is, the peak in Imx(g, uj + i0 + ) 
versus frequency is a strong function of momentum. In fact, the momentum dependence 
is so strong that it would imply the system to be very close to a magnetic instability, in 
contradiction to experiment. This problem is discussed further by Monthoux and Pines.0 
By contrast, in the RULN model the u peak position is independent of momentum. The 
latter behavior would seem unphysical from an effective RPA viewpoint, but is claimed to be 
consistent with presently available neutron scattering data on YBC0x3 On the other hand, 
it seems quite clear that substantial high frequency spectral weight is present in the undoped 
insulating cuprate parent compounds, since their AF magnon spectra must extend up to 
several hundred meV, assuming accepted estimates of the in-plane AF exchange constant 
J ~ lOOmeV. It is therefore quite plausible that additional high-frequency spectral weight 
may also exist in heretofore unexplored energy regimes in the doped cuprates. 

We thus feel that the crux of the matter lies in the neutron scattering data. The MMP 
form is based on NMR data which is essentially a zero- frequency measurement. We would 
argue that it is somewhat dangerous to infer the full frequency dependence of the suscepti- 
bility based on that data. There is hope that sufficiently reliable neutron scattering data, 
elucidating the complete momentum and frequency dependence, will become available to 
resolve these issues. We emphasize this since the significance and interpretation of current 
neutron scattering data has been a very controversial issue, especially in YBCO. 

In this context the most recent neutron scattering data by Yamada et a/.El and by Hayden 
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et alM on La2- x Sr x Cu04 (LSCO, with doping x of 14 — 15%) are of considerable interest. 
These new data seem to implyil that in the doped samples, the high-energy spectral weight 
(up to ~ 300meV) in Im% is substantially suppressed relative to the undoped, insulating 
antiferromagnetic parent (x = 0) compound and that substantial low-energy spectral weight 
(~ 22meV) is being built up by doping. Also, the width in momentum space is relatively 
independent of frequency and implies a correlation length of order a. This observation would 
be quite consistent with an RULN picture. On the other hand, low frequency data@ show 
incommensurate peaks with a correlation length of 6.7 a. This rapid change of width in 
momentum space from low to high frequencies would be consistent with an MMP picture. 
It thus appears at this point that the new LSCO data are, in some sense, intermediate 
between the RULN and the MMP models. As Hayden et al. point out,^l their data are 
consistent with a broadened magnon dispersion, the functional form of which does not look 
like either model considered here. 

On the theoretical side, we caution against over-interpreting the quantitative significance 
of results obtained from phenomenological spin fluctuation exchange models of the type 
discussed here. In order to produce a low resistivity and a high T c , the typically required 
magnitude of the Migdal parameter s in these models is sufficiently large to render the whole 
Migdal-Eliashberg approach invalid or, at the very least, makes it a rather questionable 
foundation for a quantitative theory. We also emphasize the often overlooked point that the 
conditions of applicability for the Fermi-surface-restricted Boltzmann transport theory are 
essentially the same as the conditions of validity for the Migdal approximation. Thus, as 
far as different transport formalisms are concerned, the Fermi-surface-restricted Boltzmann 
formalism is in principle "no worse" than the full- zone Bethe-Salpeter theory and, along 
with the Migdal approximation, either approach becomes quantitatively unreliable when 
the Migdal parameter s is as large as we have found it for the spin fluctuation models 
considered here. 

Finally, the required, rather large boson energy scales, with f2/e^~0.1 — 0.2, imply that 
extended (and, especially, nearest neighbor) Coulomb repulsions will substantially suppress 
the (i-wave T c in the current AF spin fluctuation exchange models, since the conventional 
mechanism of reducing the effective Coulomb interaction strength via retardationl! becomes 
largely inoperative. Calculations based on current spin fluctuation models have so far ne- 
glected extended Coulomb interactions entirely and may thus seriously overestimate the 
maximum achievable T r .0 



ACKNOWLEDGMENTS 

The authors would like to thank P. B. Allen, P. Monthoux and R. J. Radtke for several 
helpful discussions. This work was supported by the National Science Foundation (DMR 
91-20000) through the Science and Technology Center for Superconductivity (H.B.S.) and by 
the U. S. Department of Energy, Basic Energy Sciences, under Contract $W-31-109-ENG- 
38 (M.R.N.). H.B.S. would also like to acknowledge support from the National Science 
Foundation under grant No. DMR-9215123 and computing support from the University of 
Georgia and from the National Center for Supercomputing Applications at the University 
of Illinois. M.R.N, acknowledges support from the Aspen Center for Physics where part of 



19 



this manuscript was completed. 



20 



REFERENCES 



1 B. Goss Levi, Phys. Today 46, No. 5, 17 (1993). 
2 D. J. Scalapino, Physics Reports 250, 329 (1995). 

3 D. J. Scalapino, E. Y. Loh, Jr., and J. E. Hirsch, Phys. Rev. B 34, 8190 (1986); ibid. 35, 
6694 (1987). 

4 K. Miyake, S. Schmitt-Rink, and C. M. Varma, Phys. Rev. B 34, 6554 (1986). 

5 N. E. Bickers, D. J. Scalapino, and R. T. Scalettar, Intl. J. Mod. Phys. B 1, 687 (1987); 
N. E. Bickers, D. J. Scalapino, and S. R. White, Phys. Rev. Lett. 62, 961 (1989); Phys. 
Rev. B 43, 8044 (1991); C.-H. Pao and N. E. Bickers, Phys. Rev. Lett. 72, 1870 (1994); 
Phys. Rev. B 49, 1586 (1994). 

6 T. Moriya, Y. Takahashi, and K. Ueda, J. Phys. Soc. Jap. 59, 2905 (1990); K. Ueda, T. 
Moriya, and Y. Takahashi, in Electronic Properties and Mechanisms of High-T c Supercon- 
ductors, eds. T. Oguchi et ai, pp. 145 (North- Holland, 1992). 

7 A. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 167 (1990). 

8 P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. Lett. 67, 3448 (1991); Phys. Rev. 
B 46, 14803 (1992). 

9 P. Monthoux and D. Pines, Phys. Rev. Lett. 69, 961 (1992). 
°P. Monthoux and D. Pines, Phys. Rev. B 47, 6069 (1993). 
1 P. Monthoux and D. Pines, Phys. Rev. B 49, 4261 (1994). 

2 R. J. Radtke, S. Ullah, K. Levin, and M. R. Norman, Phys. Rev. B 46, 11975 (1992). 
3 Q. Si, J. P. Lu, and K. Levin, Phys Rev. B 45, 4930 (1992); Y. Zha, Q. Si, and K. Levin, 

Physica C 201, 289 (1992). 
4 R. J. Radtke, K. Levin, H.-B. Schuttler, and M. R. Norman, Phys. Rev. B 48, 15957 

(1993). 

5 J. W. Serene and D. W. Hess, Phys. Rev. B 44, 3391 (1991). 

6 R. J. Radtke, K. Levin, H.-B. Schuttler, and M. R. Norman, Phys. Rev. B 48, 653 (1993). 

7 By including extended Coulomb repulsion terms into our full-zone Eliashberg equations for 
the MMP model, Refs. [7-11], we have shown that even a small nearest neighbor Coulomb 
repulsion, of strength V 1 = 200meV, is sufficient to suppress the ci-wave T c by about a 
factor of 2. 

8 A. A. Abrikosov and L. P. Gor'kov, Zh. Eksp. Teor. Fiz. 39, 1781 (1960) [JETP 12, 1243 
(1961)]. 

9 A. J. Millis, S. Sachdev, and C. M. Varma, Phys. Rev. B 37, 4975 (1988). 

;0 W. McMillan, Phys. Rev. 167 331, (1968); P. B. Allen and R. C. Dynes, Phys. Rev. B 
12, 905 (1975). 

!1 For a comprehensive review, see P. B. Allen and B. Mitrovic, in Solid State Physics, eds. 

H. Ehrenreich et ai, Vol 37, pp. 1-91 (Academic Press, New York, 1982). 
!2 The authors would like to thank P. Monthoux for clarifying discussions concerning this 

point. 

13 For a review of LDA bandstructure calculations, see W. E. Pickett, Rev. Mod. Phys. 61, 
433 (1989). 

;4 J. C. Campuzano, G. Jennings, M. Faiz, L. Beaulaigue, B. W. Veal, J. Z. Liu, A. P. 
Paulikas, K. Vandervoort, H. Claus, R. S. List, A. J. Arko, and R. J. Bartlett, Phys. Rev. 
Lett. 64, 2308 (1990). 



21 



25 D. M. Ginsberg, W. Lee, and S. Stupp, Phys. Rev. B 47, 12167 (1993). 

26 J. Orenstein, G. Thomas, A. J. Millis, S. L. Cooper, D. H. Rapkine, T. Timusk, L. F. 
Schneemeyer, and J. V. Waszczak, Phys. Rev. B 42, 6342 (1990); M. Suzuki, in Strong 
Correlation and Superconductivity, eds. H. Fukuyama, S. Maekawa, and A. Malozemoff 
(Springer, Berlin, 1989). 

27 Both the experimental data and calculated model results for cr(iv) may already deviate 
noticeably from a simple Drude law, Eq.(7), for photon frequencies u ~ 1/ Tp. To get at least 
a meaningful order-of-magnitude estimate for a "Drude" plasma frequency uj Pi d, we have 
thus parametrized the zero-frequency peak in a{uj) in terms of its peak height cx(0) = l/p 
and in terms of its peak width, 1/td, defined here as the half-width at half-maximum 
(HWHM). We then use Eq.(13) as the defining relation for oo PtD , that is oj PtD = (An/pTo) 1 ^ 2 . 

28 P. B. Allen, Phys. Rev. B 3, 305 (1971); ibid. 17, 3725 (1978); Comments Cond. Matter 
Physics 15, 327 (1992); P. B. Allen, W. E. Pickett, and H. Krakauer, Phys. Rev. B 37, 
7482 (1988). 

29 R. Kubo, J. Phys. Soc. Japan 12, 570 (1957). 

30 In their original paper, RULN, Ref.[12], stated their coupling strength in terms of a 
dimensionless coupling parameter J 2 which, in our present notation, is given by J 2 = 
g 2 NpC/7i where Np is the non-interacting density of states at the Fermi level. For a 
hole concentration n/j = 1.18 holes per planar Cw-site (i.e. a hole doping concentration of 
x = 18%) and t = t 2 /t 1 = 0.45, we have N F = 0.2481/ti and hence a conversion factor 
N F C/ir = 0.01500eV~ 1 /ti. Thus, for an assumed t 1 = 80meV, the original upper limit 
value of I 2 = 0.10, proposed by RULN, corresponds to a value of g 2 = 0.533eV 2 . 

31 In RULN, Ref. [12], the transport relaxation time results where actually expressed in terms 
of the dimensionless transport A parameter, \ tr0 = (27rTr ) _1 , from Eq.(10), which turns 
out to be approximately T-independent in the RULN model for T~200K. RULN adjusted 
their coupling strength g 2 (or P, see previous footnote) so as to match the calculated A tf , )0 
to the measured, approximately T-independent \ tr ,D = (^Ttd)^ 1 . This procedure is of 
course entirely equivalent to adjusting g 2 so that r becomes equal to td- 

32 T. D. Holstein, Ann. Phys. (N. Y.) 29, 410 (1964). 

33 The quasi-particle mass enhancement due to spin fluctuations, Z s f, could be reduced, 
from the T = estimate Z s /~1 + \z, by thermal effects. In any boson exchange model of 
the types considered here, (including e.g. electron-phonon and spin fluctuation exchange 
systems), one can show that the relevant boson-exchange-generated mass enhancement 
factor Z (i.e. Z s f in the present context) approaches unity when the temperature T be- 
comes comparable to or larger than the characteristic boson energy scale Q [c.f., for ex- 
ample J. Zhong and H.-B. Schiittler (University of Georgia preprint, 1994)]. For the MMP 
model, Q, as estimated e.g. by the low-order frequency moments of the spin fluctuation 
spectrum a z F(u), Eq.(23) (see Table I), is of the order of 1000K and thus substantially 
larger than the relevant T range T ~ 300K. However, for the RULN model, Q is only of 
order 300 — 400K and thermal effects could thus be noticeable in the T-range of interest 
here. 

34 H.-B. Schiittler and C.-H. Pao, Phys. Rev. Lett, (in press). 

35 H. J. Vidberg and J. W. Serene, J. Low Temp. Phys. 29, 177 (1977). 

36 P. B. Allen, Phys. Rev. B 34, 4331 (1986). 



22 



A. A. Abrikosov, L. P. Gor'kov, and I. E. Dzyaloshinsky, Methods of Quantum Field 
Theory in Statistical Mechanics (Dover, New York, 1975); A. L. Fetter and J. D. Walecka, 
Quantum Theory of Many Particle Systems (McGraw-Hill, New York, 1971); G. D. Mahan, 
Many-Particle Physics, 1st ed. (Plenum Press, New York, 1981). 
J. Hertz and K. Levin, Solid State Comm. 18, 803 (1976). 
A. J. Millis, Phys. Rev. B 45, 13047 (1992). 

N. Bulut, D. J. Scalapino, and S. R. White, Phys. Rev. B 47, 2742 (1993). 
J. R. Schrieffer, J. Low Temp. Phys. 99, 397 (1995). 

K. Yamada, Y. Endoh, C.-H. Lee, S. Wakimoto, M. Arai, K. Ubukata, M. Fujita, S. 
Hosoya, and S. M. Bennington, J. Phys. Soc. Japan 64, 2742 (1995). 
S. M. Hayden, G. Aeppli, H. A. Mook, T. G. Perring, T. E. Mason, S.-W. Cheong, and 
Z. Fisk, preprint (1995). 

T. E. Mason, G. Aeppli, and H. A. Mook, Phys. Rev. Lett. 68, 1414 (1992). 



23 



TABLES 



TABLE I. Comparison of the RULN model, Eq.(4), and the MMP model, Eqs.(l) and 
(3). Abbreviations FK, VB, and BS indicate transport results obtained at T = 200K by the 
Fermi-surface-restricted Kubo and variational Boltzmann formalism, Eq.(lO), and by solution of 
the full-zone Bethe-Salpeter equation, from Ref. [11], respectively. 



Model: 




RULN 


RULN 


RULN 


RULN 


MMP 


Column: 




1 


2 


3 


4 


5 


Model input parameters: 


h [meV] 




80 


40 


250 


250 


250 


9 2 [eV 2 ] 




0.533 


0.267 


2.53 


6.82 


0.410 


/ 




1.00 


1.00 


1.00 


2.70 


1.00 


N L 




1 


2 


2 


2 


2 


n h 




1.18 


1.18 


1.18 


1.18 


1.25 


d-Wave T c and transport: 


T C [K] 




10.6 


11.2 


18.4 


99.5 


99.5 


p [fiflcm] 


FK 


447 


447 


109 


109 


109 




VB 


556 


556 


136 


150 


132 




BS 










64 


u p ,o [eV] 




0.97 


0.97 


2.42 


2.42 


2.30 


w p ,d [eV] 


BS 










1.2 


l/r [meV] 


FK 


57 


57 


86 


86 


78 




VB 


71 


71 


107 


119 


94 




BS 










46 


1/td [meV] 


BS 










12 




FK 


0.52 


0.52 


0.79 


0.79 


0.72 




VB 


0.65 


0.65 


0.99 


1.09 


0.87 




BS 










0.42 


Eliashberg A parameters: 




at T c 


1.07 


1.07 


1.64 


4.11 


1.78 




at 200K 


0.83 


0.83 


1.27 


3.41 


1.62 




at T c 


0.62 


0.62 


0.95 


2.37 


2.03 


K 


at T c 


0.30 


0.30 


0.36 


0.46 


0.73 


Frequency moments: 


(uj)z [meV] 


at T c 


32.5 


32.5 


32.8 


98 


98 




at 200K 


40.3 


40.3 


40.3 


109 


105 


(c 2 )f [meV] 


at T c 


36.9 


36.9 


37.4 


115 


141 




at 200K 


47.4 


47.4 


47.4 


128 


146 


(uj) d [meV] 


at T c 


32.5 


32.5 


32.8 


98 


71 


{" 2 )T [meV] 


at T c 


36.9 


36.9 


37.4 


115 


109 
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FIGURES 



FIG. 1. d-wave superconducting transition temperature T c versus scaled coupling strength 
g 2 jt\ for the RULN model interaction potential V(q,uj) = g 2 XRULN(l,^), Eqs.(4) and (5), with 
t = t2/t\ =0.45 and different values of t±, as indicated in the figure. T c was calculated by full-zone 
solutions of the linearized Eliashberg equations on a fc-mesh corresponding to 64 x 64 /c-points 
covering the full 1st Brillouin zone, with fermion Matsubara frequencies v m = (2m — \)irT up to 
to = 1024, using FFT techniques, as described in Refs. [9-11,14-15]. 

FIG. 2. Resistivity p versus temperature T for the MMP model, obtained from the variational 
Boltzmann (VB, dot-dashed line) and from the Fermi-surface-restricted Kubo formalisms (FK, 
dashed line) Eqs.(8-12). For comparison, we also show the results of Monthoux and Pines, MP-II, 
Ref.fll], obtained from full-zone solutions of the Bethe-Salpeter transport equation (BS, full line) 
for identical model parameters, as given in the last column of Table I, with T-dependent uj s f(T), 
X Q (T) and£(T) given by Eq. (3). 

FIG. 3. (a) ot 2 tr F(u), Eq.(ll), (b) oi 2 z F{u), Eq.(23), and (c) ct 2 d F(u), Eq.(28), versus frequency 
co for the MMP model, Eq.(l), and for the RULN model, Eq.(4), at different temperatures T. The 
RULN model parameters used are those from column 3 of Table I. The MMP model parameters 
used are those from the last column of Table I with T-dependent co s f(T), xq(T) and (,(T) given 
by Eq.(3). The 3 different temperatures include the d-wave superconducting T c values of the 
RULN model, T C:RULN = 18. 4K, and of the MMP model, T C:M mp = 99. 5K, for the respective model 
parameters, as given in Table I, columns 3 and 5, respectively. 
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